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Abstract. Recent developments in using electromagnetic meson production 
reactions to study the structure of nucleon resonances are reviewed. Possible 
future works are discussed. 



1. Introduction 

The study of nucleon resonances (N*) has long been recognized as an essential step 
towards developing a fundamental understanding of strong interactions. Since the 
discovery of the A(1232) in 1951, many iV*'s have been identified mainly through 
amplitude analyses of pion-nucleon elastic scattering. This effort peaked in the 
1970's with fairly consistent results from several independent groups. A long- 
standing challenge in hadronic physics is to understand the spectroscopy of these 
resonant states, including their electromagnetic and strong decays, within a framework 
consistent with QCD. However, many of the assumptions used in phenomenological 
analyses remain to be justified theoretically. Also, the limitations of the data base 
lead to rather large errors in many of the determined resonance parameters. 

Some additional but limited pion-nucleon scattering data have been obtained since 
the 1980's. While these data have improved the determination of ttN amplitudes at 
low energies, there has not been much impact in the N* region. The situation changed 
drastically in the 1990's with the construction of new electron and photon beam 
facilities. Experiments at Jefferson Laboratory MIT-Bates, LEGS, Mainz, Bonn, 
GRAAL, and Spring-8 are providing new data on the photo- and electroproduction 
of 7T, r), K, to, <t> mesons and 2tt final states. These data offer a new opportunity to 
understand in detail how N* properties emerge from the non-perturbative aspects of 
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QCD. For example, the extraction of 7iV — > N* transition form factors could shed light 
on the dynamical origins of the confinement of constituent quarks and the associated 
meson cloud. 

To make progress, we need to perform amplitude analyses of the data to extract 
the N* parameters. More importantly, we need to develop dynamical reaction models 
to interpret the extracted N* parameters in terms of QCD. At the present time, the 
achievable goal is to test the predictions from various QCD-based hadron models such 
as the well-developed constituent quark modelpQ and the covariant model based on 
Dyson-Schwinger Equations 2J. In the near future, we hope to understand the N* 
parameters in terms of Lattice QCD (LQCD). 

In the A region, both the amplitude analyses and reaction models have been well 
developed. We find that these two efforts are complementary. For example, a recent 
review of amplitude analyses [3] gave the result that the N-A Ml transition strength 
is Gm(0) = 3.18 ± 0.04. This value is about 40 % larger than the constituent quark 

model predictions G M (0) = |y / ^G p (0) ~ 2.30. Using dynamical models [1(5] 

one can show that the discrepancy is due to the pion cloud which is not included in 
the constituent quark model calculations. The effect of the meson cloud on the N-A 
transition has been further revealed in the study of pion electroproduction reactions. 

In the higher mass N* region, the situation is much more complicated because 
of many open channels. Any reliable analysis of the meson production data in this 
higher energy region must be based on a coupled-channel approach. In recent years, 
significant progress has been made in this direction. 

In this article, we will briefly summarize the status of theoretical and experimental 
efforts in the study of N* resonances. In section 2, we review most of the current meson 
production reaction models used to extract N* information from meson production 
data. In section 3, we summarize the results from recent analyses of photo- and 
electro-production data. The focus will be on the N — A (1232) transition form 
factors which have now been investigated in great detail. We will also summarize the 
N-N* transition form factors extracted recently by several groups. Discussions on 
future developments will be given in section 4. 



2. Theoretical developments 

Most of the models for meson production reactions can be derived by considering the 
following coupled-channel equations 

T a ,b(E) = V a , b + V a , c9c {E)T c , b {E) , (1) 

c 

where a,b,c — 7 TV, ttN, r)N, luN, KY, it A, pN aN, ■■. The interaction term is defined 
by V a . b =< a\V\b > with 

V = v bg + v R . (2) 

Here v bg represents the non- resonant (background) mechanisms such as the tree 
diagrams illustrated in Figs. DJa)-(d), and v R describes the N* excitation shown in 
Figure [lje). Schematically, the resonant term can be written as 

A*) = E^> (3) 

N* 1 
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Figure 1. Tree-diagram mechanisms of meson- baryon interactions. 

where I\ defines the decay of the z-th TV* state into meson-baryon states, and Mf is 
a mass parameter related to the resonance position. 
The meson-baryon propagator in JT]) is 

g c (E) =< c | g(E) | c > , 

with 

9 ™ = E—hve 

= g p (E)~nrS(E-H ), (4) 
where Hq is the free Hamiltonian and 

fW-E^Ho- (5) 



Here P denotes taking the principal-value part of any integration over the propagator. 
If g(E) is replaced by g p (E) and T a ^(E) by K a ^(E), |T]) then defines the K-matrix 
which is related to T-matrix by 

T a>b (E) = K a . b (E) -Y,T a ,c( E )[^ S ( E ~ H )] c K c , b (E) . (6) 

C 

By using the two potential formulation, one can cast |T]) into 

T a AE)=t b a ° b (E)+t* b (E) (7) 

with 

ta,b(E)= £ ^ N ,JE)[G(E)] hJ f N; , b (E). (8) 

NT ,Nf 

The first term of (JT]) is determined only by the non-resonant interaction 

t%{E)=v% + Y,v^{E)t%{E). (9) 

c 

The resonant amplitude ([8]) is determined by the dressed vertex 

f N*,a(E) = Tpf <a + £ rjv*,bffb (E)t b b 9 a (E), (10) 
b 

and the dressed propagator 

[G{E)-\ 3 {E) = (E- M?)S id - £^(25) . (11) 
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Here the mass-shift term is 

= E r ka^) f iV;,a(£). (12) 
a 

Note that the meson- baryon propagator g a {E) for channels including an unstable 
particle, such as 7rA, pN and aN, must be modified to include a width. 

With the above equations, we now can derive in the next few subsections most of 
the recent models for analyzing the data of meson production reactions. 



2.1. Unitary Isobar Model 

In recent years, the most widely used model is the Unitary Isobar Model (UIM) 
developed [6] by the Mainz group. This model, called MAID, is based on the on- 
shell relation ©. By including only one hadronic channel, irN (or r)N ), $Q leads 
to 

Tttn,~/N — e 1 *** cosS^nKkn^n , (13) 

where S^n is the pion- nucleon scattering phase shift. By further assuming that 
K V = v hg + v R , (|13p can be cast into the following form 

lW, 7A r( UIM) = e 5 ^cos6^ N [v% tjN ] + ]T T^ N (E) . (14) 

NT 

The non-resonant term v ha in (I14p is calculated from the standard Born terms 
but with an energy-dependent mixture of pseudo-vector (PV) and pseudo-scalar 
(PS) ttNN coupling. For the resonant term in (fT4"]) , MAID uses the following 
parameterization by Walker [7] 

r Ale^ 

T v N,jn( E ) = f*N( E ) M 2 _ ff 2 iMi T tot ^ N ^ A% ' 

where f^ N (E) and f^ N (E) are the form factors describing the decays of N* , T tot is 
the total decay width, and A 1 is the 7^V — > N* excitation strength. The phase $ is 
required by the unitarity condition and is determined by an assumption that relates 
the phase of the total photoproduction amplitude to the wN scattering phase shift. 

The UIM developed by the JLab- Yerevan collaboration [8l [9l [10] is similar to 
MAID. The main difference is that this model uses the Regge parameterization 
to define the amplitudes at high energies. Schematically, they write non-resonant 
amplitude [Background] — [e S7rN cosS 7r ]yv b9 } of (fT4")l as 

[Background] — [Born + p + lu]uim at s < Sq (16) 

1 

= [Born + p + io\uiM-——} ^7 

1 + (s - s y 

(s s )^ 

+ [Regge — Poles] ; at s > so (17) 

1 + [s - s ) 

where so ~ 1-2 GeV was determined phenomenologically in fitting the pion 
photoproduction data. 

Both MAID and JLab-Yeveran UIM have been applied extensively to analyze the 
data of 7r and r\ production reactions. 
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2.2. K-matrix Coupled- Channel Models 

2.2.1. VPI-GWU Model The VPI-GWU [TT] (SAID) can be derived from © by 
considering three channels: jN, ttN, and 7rA. The solution of the resulting 3x3 
matrix equation leads to 

TyN,nN{SAID) = Ai(l + iT n N^N) + ArT^N^N 7 (18) 

where 

At = J^jN,wN 7p , (ly) 

-"■7riV,7rA 

A R = . (20) 

The irN amplitude T^n^n needed to evaluate (|T8|) is available in SAID. In actual 
analyses, they simply parameterize Ai and Ar as polynomial functions of the on-shell 
momenta for pion and photon with adjustable parameters. Once the parameters of 
A j and Ar have been determined by fitting the pion photoproduction data, the N* 
parameters are extracted by fitting the resulting amplitude T 7 jv,7rAr at energies near 
the resonance position to a Breit-Wigner parameterization (similar to (fT5|) ). 

An extensive database of pion photoproduction experiments has been analyzed 
by SAID. The extension of SAID to the analysis of pion electroproduction data is 
being pursued. 

2.2.2. Giessen and KVI Models The K-matrix coupled-channel models developed by 
the Giessen group [H 031 HH [15] and the KVI group [16] can be obtained from §6§ 
by taking the approximation K — V . This leads to a matrix equation involving only 
the on-shell matrix elements of V 

T a , b (E) - + iViE^-X^biE) . (21) 

c 

The interaction V = v hg +v R is calculated from tree diagrams such as those illustrated 
in Figure [U The Giessen group has performed analyses [H [13l HU [15] with 7JV, irN, 
2irN, r)N, ujN, and KA channels. They further simplify 2irN as a stable partice 
channel. The KVI group [TB] has focused on the hyperon production reactions by 
performing analyses with jN, ttN, r/N, KA, KT,, cf>N channels. 

2.3. Dynamical Models 

The dynamical models of meson-baryon reactions account for the off-shell scattering 
dynamics through the use of integral equations such as ([T]) or their equivalence 
(|7l)- (fT2]) . The off-shell dynamics are closely related to the meson-baryon scattering 
wavefunctions in the short-range region where we want to map out the structure 
of N*. Thus the development of dynamical models is an important step towards 
interpreting the extracted N* parameters. 

In recent years, the predictions from the dynamical models of Sato and Lee (SL) 0] 
and the Dubna-Mainz- Taiwan (DMT) collaboration [5] are most often used to analyze 
the data in the A region. The SL model can be derived from (f7])- ([T2"T) by keeping 
only one resonance N* — A and two channels a,b = ttN, 7V. In solving exactly the 
resulting equations the non-resonant interactions v^ 9 N vN and v^ N lN are derived from 
the standard PV Born terms and p and to exchanges by using an unitary transformation 
method. 
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The DMT model also includes only irN and 7 A channels. They however 
depart from the exact formulation based on ((T|) or ([?]) - ([!!?]) by using the Walker 
parameterization (|15p to describe the resonant amplitude. Accordingly, their definition 
of the non-resonant amplitude also differs from that defined by ([9]) : t b c 9 b in the right- 
hand side of §§§ is replaced by the full amplitude T c ^. Furthermore, they follow MAID 
to calculate the non-resonant interaction v^f N jN from an energy-dependent mixture 
of PS and PV Born terms. 

Extensive data of pion photoproduction and electroproduction in the A region 
can be described by both the SL and DMT models. However, they have significant 
differences in the extracted electric E2 and Coulomb C2 form factors of the -fN — > A 
transition. Both models show very large pion cloud effects on the 7 A — > A transition 
form factors in the low Q 2 region. 

The Utrecht-Ohio model [17] has also succeeded in describing the data in the A 
region. Despite some differences in treating the gauge invariance problem and the A 
excitation amplitude, its dynamical content is similar to that of SL and DMT models. 

The equations ([7|- (|12p are used in a 2- A* and 3-channels (ttN, rjN, and irA) 
study [18] of irN scattering in Sn partial wave. This work illustrated the extent to 
which the quark-quark interactions in the constituent quark model can be determined 
directly by the ttN reaction data. The equations ([7|)- (fT2]) have also been used to 
show that the coupled-channel effects due to the ttN channel are very large in u> 
photoproduction [TH] and K photoproduction [5D]. All A* identified by the Particle 
Data Group are included in these two dynamical calculations. 

The coupled-channel study of both irN scattering and jN — > irN in the Sn 
channel by Chen et al [21] includes irN, r]N, and 7A channels. Their ttN scattering 
calculation is performed by using ([1]). For their 7 A — > nN calculation, they neglect 
the 7./V — > nN — > irN channel coupling, and follow the procedure of the DMT model 
to evaluate the resonant term in terms of the Walker parameterization (| 1 5|) . They 
find that four A* are needed to fit the empirical amplitudes in the Sn channel up to 
W = 2 GeV. This approach is being extended to also fit the nN — ► nN and 7 A — > nN 
amplitudes in all partial waves. 

A coupled-channel calculation based on ([T]) has been carried out by the Julich 
group [55] for nN scattering. They are able to describe the irN phase shifts up 
to W = 1.9 GeV by including the irN, rjN, irA, pN and crN channels and 5 N* 
resonances: P 3 3(1232), S u (1535), S u (1530), 5i 3 (1650) and L>i 3 (1520). They find that 
the Roper resonance Pi 1 (1440) arises completely from the meson-exchange coupled- 
channel effects. 

A coupled channel calculation based on ([T]) for both the irN scattering and 
7 A — > ttN up to W = 1.5 GeV has been reported by Fuda and Alarbi [23]. They 
include the nN, 7A, -qN and ttA channels and 4 N* resonances: P 33 (1232), Pn(1440), 
5*11 (1535) and Di 3 (1520). They adjust the parameters of their model to fit the 
empirical multipole amplitudes in low partial waves. 

More simple coupled-channel calculations have been performed by using separable 
interactions. In the model of Gross and Surya [24] . separable interactions come from 
simplifying the meson-exchange mechanisms in Figs [IJa)-[Hc) as a contact term like 
Figure [Hd). They include only irN and 7A channels and 3 resonances: P 33 (1232), 
Pn(1440) and L>i 3 (1520), and restrict their investigation up to W < 1.5 GeV. To 
account for the inelasticities in Pn and D13, the N* — > irA coupling is introduced in 
these two partial waves. The inelasticities in other partial waves are neglected. 
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A similar separable simplification is also used in the chiral coupled-channel 
models [U] [53] for strange particle production. There the separable interactions are 
directly deduced from the SU(3) effective chiral lagrangians. They are able to fit 
the total cross section data for various strange particle production reaction channels 
without introducing TV* resonance states. In recent years, this model has been further 
extended by Lutz and Kolomeitsev [27] to also fit the irN scattering data. 

None of the dynamical models described above account for all of the effects 
arising from the irirN channels, which contribute about 1/2 of the irN and 7 A 
total cross sections in the higher mass N* region. Consequently, the N* parameters 
extracted from these models could have uncertainties due to the violation of the ttttN 
unitarity condition. One straightforward way to improve the situation is to extend the 
Hamiltonian formulation of (4] to include (a) p — » tttt and a —> tttt decay mechanisms 
as specified in Ty of ([24]) . (b) for non-resonant tttt interactions, (c) vmn,-kitN for 
non-resonant MN — > ttttN transitions with MN = jN or ttN, and (d) v^-kN^-tcN f° r 
non-resonant ttttN — > ttttN interactions. Such a dynamical coupled-channel model 
has been developed recently by Matsuyama, Sato, and Lee [28] (MSL). The coupled- 
channel equations from this model can also be cast into the form of (j7j)- (fl"2|) except 
that the driving term of ([9]) is replaced by 

< 3 , b ^V a , b = vl% + Z a , b (E), (22) 

where 

Z a>b (E) = < a I F— — ^ —Ft I b > 

hi — Ho — V n N,TrN — Vtttt — VkkN^ttN + «e 

- <a\Z v (E)\b> S a , b (23) 

with 

= Ta^ttW + hp t7I7T + h a ,TT-K + V„N,1T-kN ■ (24) 

Note that Sy(F) in ([23]) describes the self-energy of unstable particle channels 
7rA, pN, and aN . Obviously Z a ^[Fi) contains the non-resonant multiple scattering 
mechanisms within the ttttN subspace. It generates ttttN unitarity cuts. It has been 
shown in [28 that this term has very large effects on 2n production reactions. This 
dynamical coupled-channel model is currently being further developed at the Excited 
Baryon Analysis Center (EBAC) at JLab. 

3. Experimental developments 

The past ten years has seen an unprecedented increase in the volume and precision of 
N* data from electromagnetic meson production reactions. This has resulted from 
the development of CW electron accelerators, which drastically reduce accidental 
backgrounds in coincidence experiments, as well as from advances in charged, neutron 
and photon detector technologies such as nearly An solid angle coverage and improved 
spatial/energy resolution. In addition, the development of polarized beams, targets 
and polarimeters has opened up a previously unavailable set of observables which 
are sensitive to the interference between resonant and non-resonant processes in the 
scattering amplitude. An extensive survey of initial data from the new facilities can 
be found in the review by Burkert and Lee [31] . In this section we highlight results 
from recent analyses of ongoing experiments. 
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Figure 2. W dependence of partial wave coefficients from Legendre fit of 
experimental n° electroproduction structure functions <ry_|_£,, u^t an d &tl- Solid 
curve shows UIM fit results. Vertical line shows position of A(1232) resonance. 
CLAS data [33] taken at E=f.046 GeV and Q 2 = 0.16 GeV 2 . 

3.1. Resonance transition form factors from it, rj, 2tt production data 

Pion electroproduction provides new information not available using real photons. 
First, measuring the Q 2 dependence of the photocoupling amplitudes probes a 
range of distance scales within the nucleon, revealing spatial information about 
the constituents. Second, polarization of the virtual photon allows study of both 
longitudinal and transverse modes of resonance excitation. Both of these factors place 
strong constraints on the structure of the transition current used in phenomenological, 
quark model and lattice QCD calculations. 

3.1.1. A(1232) Resonance A longstanding question is the origin of the non-zero 
quadrupole strength experimentally observed in the j*p — > A(1232) transition. 
Quadrupole transitions, which occur via the absorption of E2 or C2 photons, imply 
non-zero orbital angular momentum in the N or A wave function, which can arise from 
configuration mixing due to residual gluon or meson interactions between constituent 
quarks or from the direct photon coupling to the pion cloud. Until recently the data 
were not sufficiently accurate to constrain calculations. New precision measurements 
of the p(e, e'p)ir° reaction at BATES, JLAB and MAMI have drastically improved our 
knowledge of the Q 2 dependence of the electric and Coulomb quadrupole/magnctic 
dipolc ratios Rem = E\ + /M\ + and Rsm = S\ + /M 1+ . The large acceptance CLAS 
experiments at JLAB measure the interference terms Re{E\+M± + ) and Re(Si+M* + ) 
via the azimuthal and partial wave decomposition of the pir° final state: 

= £ Aj+ L P e (cos 61) + sin 2 91 Cr cos 2# 
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Figure 3. Magnetic dipole transition form factor G* M for 7*JV — » A(1232), 
normalized to the elastic dipole Gd- Inclusive data are from pre-1990 experiments 
at Bonn, DESY and SLAC. Other data points were obtained from Unitary Isobar 
Model analysis of the exclusive p(e, e'p)7r° reaction. Curves from Sato-Lee model 
discussed in text. 



1 

+ sin6>;^Dj L P £ (cos6i;)cos0; (25) 

1=0 

where to first order, retaining only M\ + dominated terms 

|M 1+ | 2 =< +L /2, (26) 

Re{E 1+ M* l+ ) = (A T 2 +L - 2 C TT / 3)/8, (27) 
Re{S l+ M* l+ ) =Dj L /%. (28) 

Here "T" and "L" refer to the transverse and longitudinal components of the virtual 
photon polarization. The data are shown in Figure [2] for Q 2 = 0.2 GeV 2 where it 
is seen that all the Legendre coefficients in (1231) show a W dependence consistent 
with interference terms involving the strong resonant M\ + . At higher Q 2 , where M\ + 
dominance is no longer valid, corrections to (|26ti28|) can be made using the Unitary 
Isobar Model (UIM) described in Section 2.1 where non-resonant backgrounds and 
tails from higher resonances are parameterized using (fTi)) and (fT"5)) . A typical UIM fit 
is shown in Figure O 

Figure [3] shows the Q 2 dependence of the magnetic dipole G* M , which is related 
to the isospin / = 3/2 irN multipole Im(M^) through the A width and a 
kinematic factor. These new exclusive measurements of G* M now extend over the 
range Q 2 — 0.06 — 6.0 GeV 2 , and confirm the rapid Q 2 falloff relative to the elastic 
dipole seen in previous inclusive measurements, but with much greater sensitivity to 
the resonant Mi+ multipole at higher Q 2 . As previously noted, at least part of this 
rapid falloff may arise from the 'dressing' of the 7* TV — » A vertex by rescattering 
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Figure 4. Quadrupole -yN — > A(1232) transition form factors at low Q 2 (left) 
and high Q 2 (right) from recent experiments I31| . Electric quadrupole -E14- and 
scalar quadrupole 5i+ are normalized to the magnetic dipole Mi+. Preliminary 
data from CLAS [33] [35] and MAMI [40] [33] are also included. 



through the pion cloud. The 'bare' curve in Figure [3] is obtained in the SL dynamical 
model when the predicted 'dressed' G* M form factor is fitted to the data. It shows the 
increasing contribution from mesons as Q 2 — > 0. Similar results have been obtained 
in the Dubna-Mainz- Taipei model. 

The experimental quadrupole transition ratios REM and RSM are shown in 
Figure 0] The JLAB results are based on a full partial-wave analysis, while the 
MAMI and BATES data sets cover a more limited angle range. For the latter, 
MAID, DMT and SL model fits are used to estimate the A resonant multipoles. 
The other data points with errors are quenched lattice QCD calculations [30] which 
have undergone a linear chiral extrapolation to the physical pion mass. It is evident 
that while the calculations qualitatively account for the magnitude of REM and RSM, 
the Q 2 dependence is not consistently described. The low Q 2 underprediction of RSM 
by LQCD is particularly striking. It was recently demonstrated, using a relativistic 
chiral effective field theory (xEFT) calculation [32], that RSM may be particularly 
sensitive to a non-analytic dependence of pion loop diagrams on the quark mass, 
rendering invalid the linear chiral extrapolation used in [30] . A strong negative slope 
near Q 2 = is predicted by xEFT for both REM and RSM (see Figured]). While the 
lowest Q 2 RSM points from MAMI and BATES appear to agree with this prediction, 
the data for higher Q 2 are in better agreement with the SL and LQCD predictions. 
For REM, the data do not show the strong Q 2 dependence predicted by %EFT near 
Q 2 = 0, and also lie outside of the estimated theoretical uncertainty arising from 
higher-order effects. 

At the largest Q 2 measured by CLAS (right panel in Figure 2]) neither REM 
nor RSM exhibit the helicity conserving scaling behavior expected in the hard pQCD 
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Figure 5. Response functions for p(e, eff)ir° measured by Kelly etal 1441 1451 at 
W = 1.23 GeV and Q 2 = 1.0 GeV 2 . Notation refers to transverse (t), normal (n) 
or longitudinal (1) components of the proton recoil polarization. MAID, DMT, 
SAID and SL models given by dash-dotted, dotted, short-dashed and long-dashed 
curves. 



limit (REM^ 1, RSM— > constant). Normally the transverse momenta of quarks are 
expected to be of the order of Aqc d , and thus only important for resonance excitations 
at low Q 2 . However Idilbi, Ji and Ma [36] have shown that the orbital motion of 
small— x partons makes a non-negliblc contribution to helicity non-conservation in the 
NA transition, thus invoking both soft and hard QCD processes. 

In addition to understanding the long-range chiral dynamics which are evidently 
important near Q 2 = 0, there has long been a fundamental interest in the 'bare' 
jN — > A vertex, which is presumably more sensitive to the short-range interquark 
potential relevant to QCD-motivated models. The bare couplings, which also contain 
information about the quark core wave functions of the N and A, can be estimated 
in dynamical reaction models by treating them as free parameters to be fitted to the 
experimental data. This was performed using the SL model for all recent data sets 
below Q 2 = 2.0 GeV 2 . A typical fit to CLAS and MAMI data at Q 2 = 0.2 GeV 2 is 
shown in Figure while the Q 2 dependence of the extracted bare couplings g m ,g e 
and g c are shown in Figure [7] The simplest assumption, that the bare electromagnetic 
couplings have an identical Q 2 dependence, is clearly ruled out by the data. Unlike 
g m , both g e and g c show a strong departure from smooth behavior at low Q 2 . This 
behavior may indicate the presence of additional long-range physics not included in 
the reaction models. It is certain that more experimental data at the lowest Q 2 would 
aid the interpretation of these results. 

All reaction models make implicit assumptions about the mechanisms responsible 
for non- resonant backgrounds. This often leads to substantial model dependence 
in the extraction of resonance parameters due to ambiguities in relative phases, 
omitted diagrams, etc. In addition these same models are often used to perform the 
partial wave analysis of the experimental data, which should in principle be a model 
independent procedure. Clearly, more complete experiments would help resolve these 
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Figure 6. Fits of Sato-Lee dynamical model to experimental p(e, e'p)n° structure 
functions &t+l, TT and oxt- Preliminary data from CLAS (solid) |33) and 
MAMI (open) |34l are shown at Q 2 = 0.2 GeV 2 . The 'baseline' curve refers to 
the default Q 2 parameterization of the 'bare' couplings gm,9e and g c prior to the 
fit. 
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Figure 7. Q 2 dependence of bare magnetic dipole, electric quadrupole and 
coulomb quadrupole couplings g m , g e and g c for the 7*p — > A transition, extracted 
from fits of SL dynamical model to BATES, CLAS and MAMI data. Solid line 
assumes Q 2 dependence of g e and g c couplings are identical to g m . 
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ambiguities. A recent JLAB/Hall A experiment used a polarized beam and a high 
resolution magnetic spectrometer instrumented with a recoil polarimeter to measure 
16 neutral pion electroproduction response functions at Q 2 = l.O GeV 2 . Twelve of these 
response functions were observed for the first time. These data are shown in Figure[5Jat 
the peak of the A(1232), compared to MAID, DMT, SAID and SL model predictions. 
Although the overall agreement is good, there are several response functions for which 
the models clearly diverge in their predictions. This could be due simply to lack 
of previous experimental information on the weaker underlying multipoles. Due to 
the phase sensitivity of these polarization observables, a nearly model-independent 
extraction of REM and RSM was obtained. 

3.1.2. Pn(1440) A significant new result using the large acceptance data from CLAS 
is the demonstration [41] of a large sensitivity to the and multipoles. These 
multipoles receive large contributions from excitation of the 1 — 1/2 Pn(1440) or 
"Roper" resonance. Due to the isospin coupling, this sensitivity is maximized in the 
7J7T+ channel. It turns out that UIM fits to combined data sets which include the 
helicity dependent asymmetry A' LT (nn + ) and the cross section da(nn + ) taken in the 
second resonance region, significantly constrain the real and imaginary parts of these 
multipoles through their interference, based on the UIM analysis, with largely Born 
dominated non-resonant multipoles. The sensitivity of the quantity o~' LT w da ■ A' LT 
to the imaginary (resonant) part of these multipoles is demonstrated in the left panel 
of Figure [5] The plot shows the sensitivity of this observable at Q 2 = 0.4 GeV 2 
to a 0.5 [ib 1 ! 2 shift in the transverse and longitudinal resonant photocouplings. 
The sensitivity is the result of the interference term that mixes real and imaginary 
amplitudes 

a LT > ~ Im{L) ■ Re(T) - Im(T) ■ Re(L) , (29) 

where L and T represent the longitudinal and transverse amplitudes, respectively. At 
much higher Q 2 , the tail from the A(1232) is strongly diminished, and the Roper is 
directly visible as a small bump in the total nir + cross section near W — 1.4 GeV 
(right panel). 

Using UIM and MAID model fits to CLAS data, the Roper photocoupling 
amplitudes have now been extracted over a large range of Q 2 , as shown in Figure [5] A 
few of the Q 2 points have been further analyzed in a multi-channel global fit [42] which 
includes the pn and ptt + tt~ final states and the agreement with the single channel 
analysis is good. We also show recently published [131 US] points at Q 2 = 1.0 GeV 2 
from Hall A. The new CLAS06 data at even higher Q 2 [46 confirm the zero-crossing of 
A^ 2 first reported in [42] . This feature can completely rule out some previous quark 
model calculations, such as most non-relativistic approaches, as well as the quark- 
gluon hybrid approach for which no zero crossing occurs and for which S^ 2 = 0. The 
strong longitudinal response seen in these data could revive the traditional 'breathing 
mode' description of the Roper, possibly augmented by coupling to the qq cloud of 
the nucleon as suggested by the good agreement with the VMD model of Cano and 
Gonzales [49j. Dynamical models and lattice QCD calculations could provide useful 
input towards intepreting the zero-crossing seen in , 2 . 

3.1.3. Second and Third Resonance Regions The most prominent states for W > 
1.5 GeV are 511(1535), Z?i3(1520), and Fis(1680). These have fairly unambiguous 
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Figure 8. CLAS measurement of the W dependence of a LT i (left) and ctt+l 
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(right) model sensitivity to changes in multipoles M,_ and SJ_ involved in 
formation of the Roper resonance. Data were taken at Q 2 = 0.4 GeV 2 (left) and 
Q 2 = 1.72 - 4.16 GeV 2 (right). 




Q ! (GeV a ) 



Cono (71 claud) 
Copstlck (RQM) 
Simula (RQMf) 
Dorg (jiq") 




Q' (GeV ! ; 



Cono (,> Cloud) 

Cops-.ick {RQM; 

Simula [ROW) 

Dong (nq 1 ) 



Figure 9. Q 2 evolution of the Pn(1440) (Roper) amplitudes A^, 2 and Sy 2 - 
Both CLAS/UIM and Hall A points include model error. The MA05 analyses [5B] 
are dominated by pir° data. Other points are GW02 1701 . The curves are model 
calculations !49 53 55 52 [54] [50] described in the text. 
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Figure 10. Q 2 evolution of the •yp — > 5n(1535) photocouplings amplitudes -Ai / 2 
and Si/2- See Figure [9] for explantion of notation. 



signatures in partial wave analyses and can be reasonably well isolated in single pion 
production channels. Furthermore, except for the Sn, the -kN branching fractions and 
widths are known sufficiently to extract photocouplings with r* 10% accuracy. On the 
other hand, a consistent quark model description of the Q 2 dependence of these three 
orbital excitations has proven elusive. The new preliminary CLAS p(e, e'ir + )n data 
(labeled CLAS06 in Figures I51 TT2"]) has extended the study of N — > N* transition form 
factors into the unexplored region above Q 2 = 2 GeV 2 . 

For the 5n(1535), the CLAS/UIM determination [42] of photocoupling 
amplitudes from it and r\ electroproduction channels (assuming branching fractions 
(3 7:N = 0.4, = 0.55 and T to t = 150 MeV) give results which are in good agreement 
up to the highest Q 2 measured, unlike at Q 2 = where A^, 2 from it production 
is smaller than for r] production (see Figure [T0|) . It has been shown J47J [48] that 
strong channel- coupling effects complicate the analyses near Q 2 = 0, although the Q 2 
evolution of these effects has not been studied. The new data at higher Q 2 show a 
more dipole-like form factor compared to the anomalously hard form factor seen for 
Q 2 < 2 GeV 2 . Also for the first time a reliable measurement of the longitudinal 
photocoupling has been made, which may shed more light on the role of meson 
rescattering. 

A more precise determination of the 7*p — > Z?i3(1520) and 7*p — > Fis(1680) 
transitions has resulted from the increased I = 1/2 sensitivity provided by the CLAS 
7J7T+ data. The helicity non-conserving transverse amplitude A^ 2 for both transitions 

is clearly dominant at Q 2 = (see Figs. [TT] and [T2|) . and falls very rapidly with 
increasing Q 2 . All quark models predict that at high Q 2 , the helicity conserving A^ 2 
will dominate the transition. However neither the harmonic oscillator potential nor 
hypercentral potential calculations shown in the figures describe the Q 2 dependence 
satisfactorily. This is especially true for the longitudinal couplings (see also Figure [T3|) . 
where the models fail to consistently predict even the sign. For example, quark model 
wavefunctions obeying SU(6) ® 0(3) symmetry predict 5^ 2 (5n) w —S^ 2 (Di 3 ) [5T] . 
in disagreement with the data. 

Most higher lying states (M > 1.65 GeV) have large 2tt branching fractions and 
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Figure 11. Q 2 evolution of the jp —> Di3(1520) transverse photocoupling 
amplitudes Ax/2 an< ^ -^3/2 ■ See Figure [9] for explantion of notation. 
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Figure 12. Q 2 evolution of the ■yp — > _Fi5(1680) transverse photocoupling 
amplitudes A^/i and A 3/ / 2 - See Figure[9]for explantion of notation. 



reliable photocoupling amplitudes can only be obtained by performing (at minimum) 
fits to measurements of both single and double pion production channels. Also, 
detection of two pions permits the study of possible 'missing' resonances which 
have weak couplings to single meson channels. These broad, overlaping states in 
the mass region W > 1.8 GeV are predicted to decay to both Air and Np, which 
makes the 7r + 7r~ channel ideal for their detection. However, the reliable extraction 
of iV* photocouplings is also hampered by a poor understanding of the dynamical 
mechanisms behind 2tt production. Currently the most accurate data exists in the 
second resonance region, where it may be easier to disentangle the numerous resonant 
and non-resonant subprocesses which contribute to the 2tt channel. 

The analyses (SH [60j [61] of two pion production by using the JLab-Moscow State 
University (JLAB-MSU) isobar model considered only the minimum set of the tree 
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Figure 13. Q 2 evolution of the longitudinal photocoupling amplitudes S\i% f° r 
the ■yp — ► Z?i3(1520) and -yp — > Fi5(1680) transitions. See Figurel9lfor explantion 
of notation. 




Figure 14. Left: Total pw 7T electroproduction cross sections from 
CLAS [59]. Solid line: Fit using the JM05 isobar model [57]. Dashed 
line: JM05 with new 3/2 + (1720) candidate state taken out. Right: 
Q 2 dependence of photocouplings for the D33(1700) (top) and Pi3(1720) 
(bottom) states. Open circles are results extracted from CLAS 2tt data |57] . 
Other points are previous In analyses. 



diagrams proposed in the original work of Luke and Soding [55] ■ However, they made 
two improvements. They included all 3-star and 4-star resonances listed by the Particle 
Data Group and used the absorptive model to account for the initial and final state 
interactions. They found that the 7riVA form factor is needed to get agreement with 
the data of 7p — > 7r~A ++ , while the initial and final state interactions are not so large. 
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In analyzing the two-pion electro-production data, they further included a mtN phase- 
space term with its magnitude adjusted to fit the data. This term was later replaced 
by a phenomenological particle-exchange amplitude which improves significantly the 
fits to the data. Evidence for intermediate isobar channels such as 7r + Di3(1520) and 
possibly 7r + Fi 5 (1685) and 7r _ P 33 (1600) are suggested by these fits. Most prominently, 
for increasing Q 2 , the 2ir data show a peak emerging at W « 1.7 GeV (Figure [T4|) 



which is identified [BTj as a iV*(| + ,1720) state. However, this state differs from the 
conventional Pi3(1720), having a much stronger photocoupling strength and larger 
nA/pN decay ratio compared to PDG. The answer to whether this is a new state or 
a better experimental determination of an existing resonance awaits further analysis. 

Polarization can play an important role in the identification of specific 
subprocesses which contribute to 2tt production. A good example is beam-helicity- 
asymmetry data in the jp —> pn + Ti~ reaction measured by Strauch et aL|58j using 
CLAS with a circularly polarized photon beam in the mass range W = 1.35—2.30 GeV. 
The asymmetry displayed in Figs. I15ITB1 is defined as 



(30) 



where P 7 is the circular polarization of the photon and are the cross sections for the 
photon helicity states A 7 = ±1. Figure [TBI shows a Fourier analysis of the asymmetries 
measured in the mass region of the iV(1520), compared to model calculations 
incorporating ttA and pN quasi-two-body intermediate states. The A(1232) bump 
seen in the M(pir~) invariant mass (FigurdTB]) for the a% coefficient is suggestive of 



the sequential decay 2V(1520) -» tt + A c 



p. Previous unpolarized cross section 



measurements showed strong A signatures only in the M(pn + ) distributions, from 
the isospin favored D13 — > 7r~A ++ transition. This again illustrates the sensitivity of 
polarization to the relative phases between various amplitudes, which can selectively 
amplify combinations of interfering background and resonance amplitudes. 
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Figure 15. Helicity asymmetries 1581 integrated over CLAS acceptance at 
W = 1.50 GeV. Each panel indicates the value of the invariant mass M(p7r+). 
The solid, dashed curves are calculations by Mokeev et al. |6D1 [6T| and Fix and 
Arenhovel I63| . respectively. Inset shows the definitions of the azimuthal angle 4> 
for the reaction 7p — > p7r + 7r~ . 
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Figure 16. Fourier coefficients a\ and 0,2 of the helicity asymmetry as a 
function of the invariant mass M(p7r — ) for W = 1.520 GeV. The curves are the 
results of Fix and Arenhovel |63| . The diagram illustrates the sequential decay 
7V(1520) — * 7T+A — » 7r+7r _ p which might be observed in the data. 

3.2. Analyses of uN , KA, KY> production data 

Because the production cross sections of the ujN , KA, ATE channels are much weaker 
than those of the irN and 2ttN channels, tree-diagram models cannot be used reliably 
to extract the N* information from the data of photo- and electro-production of these 
weak channels. We therefore review only the results from coupled-channel analyses. 

Most extensive analyses [T4j \TE[ [16] of KA, ATE and loN photoproduction 
data have been performed by using the K-matrix coupled-channel models. The 
dynamical coupled-channel model has only been applied to analyze [20j [67] the 
KA photoproduction. The main focus of these analyses is to identify possible new 
resonances and also to extract more accurate information on N* parameters. 

The loN photoproduction has been analyzed by the Giessen group with the K- 
matrix coupled-channel model, defined by PTj) . with 7 AT, irN, r]N, 2irN, ojN channels. 
The 2wN channel is approximated as a stable particle channel with its coupling to the 
irN channel constrained by the empirical partial- wave amplitudes |69) of irN — > ttttN . 
The interaction V of (p?Tj) is defined by a set of phenomenological Lagrangians with 
the parameters and associated form factors adjusted to fit the empirical amplitudes 
of 77V — » irN, the cross sections of irN, — > and irN, — > uiN. 

In Figure[l7l we show their fits to the 7^V — > ujN differential cross sections. They 
found a strong contribution from the Dis(1675) resonance to nN — > luN and from 
the Fi5(1680) resonance to 7^ — > uiN . From their analysis in Ref [T3], they also 
extracted the — > A^* helicity amplitudes. Their results are listed in the first rows 
of Tables [TT21 The SAID values from the VPI-GWU single-channel analysis are shown 
in the second rows. The differences between these two analysis results indicate the 
importance of including the coupled-channel effects in extracting the N* parameters. 

Experimental results for KA and KT, production have become very extensive 
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over the past decade. However the coupled-channel analyses of these data are only 
now becoming available. The Giessen group, by extending their model to include KY 
channel, have found that 5n(1650), Pi3(1720) and Pi3(1900) are the main resonance 
contributions to K A photoproduction. They show that a peak near 1.9 GeV, which 
was identified with the "missing" Di3(1960) in a tree-diagram model analysis [71 j . 
can be described by a coherent sum of the contributions from the " known" resonant 
and non-resonant amplitudes. This clearly demonstrates the importance of using a 
coupled-channel approach to analyze the data of any channel, such as the KY and 
ojN channels, having much smaller production cross sections than the ttN and 2irN 
channels. 

The KVI group has analyzed the KA and KY, photoproduction data using a K- 
matrix coupled-channel model with 77V, irN, r/N, KA, KY,, and <j>N channels. They 
use the SU(3) Lagrangian to define the interaction V of (f2Tj) . With appropriate choices 
of form factors, they found that the data can be described very well with most of the 
coupling parameters not too different from their SU(3) limit. Their fits to the data 
on KA and K E total cross sections are shown in Figure HH They found very large 
contributions from both resonances and coupled-channel effects. 

The KA photoproduction has also been analyzed by using a dynamical coupled- 
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Table 1. Helicity amplitudes of / = ^ resonances ( in 10~ 3 GeV _1//2 ) determined 
in the Giessen study (first line); second line: results of SAID group 1701 ; "n.g.": 
not given. 
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Table 2. Extracted helicity amplitudes of / = | resonances ( in 10 -3 GeV -1 / 2 ) 
considered in the Giessen study. Second line: results of SAID group 1701 ; "n.g.": 
not given. 
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A l/2 
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F 35 (1905) 


32 
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2(5) 


-56(5) 
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-23 




4(6) 
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channel model [67) which is based on an extension of the dynamical model of [1] 
to include the KA channel. It differs from the K-matrix coupled- channel models of 
Giessen and KVI groups in accounting for the off-shell scattering dynamics. This 
analysis finds that the most relevant known resonances in 7p — > K + A are <Sii(1535), 
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W [GeV] 

Figure 18. Experimental data on KA and KT, total cross sections from 
SAPHIR 64, 65 compared to calculations from KVI group |16l . Dashed line 
shows effect of turning off all resonance contributions. Dotted line shows coupled 
channel effects switched off. Solid line is full calculation. 




1.75 2 2.25 1.75 2 2.25 
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Figure 19. Comparison of the W spectra of KA production from CLAS |68) 
at various CM. angles with model calculations from |67| . Solid curves: full 
calculation. Dotted curves: off-shell scattering turned off. Dashed curves: 
coupled-channels turned off. 



Pi3(1900) and Z?i3(1520). In addition, they also demonstrate the importance of off- 
shell effects and meson cloud effects in interpreting the extracted N* parameters, as 
shown in Figure [T9l 
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4. Future developments 

In the past few years, we have witnessed very significant progress in the study of N* 
physics. We now have fairly extensive data for 7T, rj, K, u>, <f> and tttt production 
channels. Many more data will soon be available. Theoretical models for interpreting 
these new data and/or extracting the N* parameters have also been developed 
accordingly. 

From the analyses of the single pion data in the A region, quantitative information 
about the 7^V — > A transition form factors have been obtained. With the development 
of dynamical reaction models, the role of pion cloud effects in determining the A 
excitation has been identified as the source of the long-standing discrepancy between 
the data and the constituent quark model predictions. Moreover, the Q 2 -dependence 
of the 77V — > A form factors has also been determined up to about Q 2 ~ 6 GeV 2 . The 
Magnetic Ml, Electric E2 and Coulomb C2 j*N — > A transition form factors should 
be considered along with the proton and neutron form factors as benchmark data for 
testing various hadron models as well as Lattice QCD calculations. 

The analyses of the it and 77 production data had led to a rather quantitative 
determination of several N* parameters in the second resonance region. However, 
a correct interpretation of the N* parameters in terms of current hadron-model 
predictions requires a rigorous investigation of the dynamical coupled-channel effects 
which are not included in amplitude analyses based on either the K-matrix isobar 
model or the dispersion relation approach. 

There has significant progress in analyzing the data for KA, ifE, and uN 
photoproduction. It is now generally agreed that these data must be analyzed in 
a coupled-channel approach. New information on N* parameters have been extracted 
from K-matrix coupled-channel analyses. The importance of both coupled- channel 
and meson cloud effects have also been revealed. It will be interesting to extend these 
analyses to analyze the electroproduction data for also extracting the N-N* transition 
form factors. 

The analysis of 2irN production data is still in a development stage, while some 
interesting results have been obtained from the analyses based on tree-diagram isobar 
models. A reliable analysis of 2wN data must satisfy the ttttN unitarity condition 
rigorously. Progress in this direction has been made recently. The dynamical coupled- 
channel model developed in Ref. [28] accounts for all relevant channels and the effects 
due to the ttttN unitarity condition. It has been applied to reveal important dynamical 
features due to the ttttN unitary cuts which must be included in extracting the N* 
information from 27r production data. The analyses based on this dynamical coupled- 
channel approach are being performed at the Excited Baryon Analysis Center (EBAC) 
at Jefferson Laboratory. 

Progress made in the past decade has resulted from close collaborations between 
experimentalists and theorists. With much more complex data to be analyzed and 
interpreted, such collaborations must be continued and extended in order to bring the 
study of N* physics to a complete success. 
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